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ABSTRACT 

Temperature  measurement  is  one  of  the  key  quantifies  in  ocean  research.  Temperature  variations  on  small  and  large 
scales  are  key  to  air-sea  interactions  and  climate  change,  and  also  regulate  circulation  patterns,  and  heat  exchange.  The 
influence  from  rapid  temperature  changes  within  microstructures  are  can  have  strong  impacts  to  optical  and  acoustical 
sensor  performance.  In  this  paper,  we  present  an  optical  fiber  sensor  for  the  high-resolution  and  high-speed  temperature 
profiling.  The  developed  sensor  consists  of  a  thin  piece  of  silicon  wafer  which  forms  a  Fabry-Perot  interferometer  (FPI) 
on  the  end  of  fiber.  Due  to  the  unique  properties  of  silicon,  such  as  large  thermal  diffusivity,  notable  thermo-optic  effects 
and  thermal  expansion  coefficients  of  silicon,  the  proposed  sensor  exhibits  excellent  sensitivity  and  fast  response  to 
temperature  variation.  The  small  mass  of  the  tiny  probe  also  contributes  to  a  fast  response  due  to  the  large  surface-to- 
volume  ratio.  The  high  reflective  index  at  infrared  wavelength  range  and  surface  flatness  of  silicon  endow  the  FPI  a 
spectrum  with  high  visibilities,  leading  to  a  superior  temperature  resolution  along  with  a  new  data  processing  method 
developed  by  us.  Experimental  results  indicate  that  the  fiber-optic  temperature  sensor  can  achieve  a  temperature 
resolution  better  than  0.00 1°C  with  a  sampling  frequency  as  high  as  2  kHz.  In  addition,  the  miniature  footprint  of  the 
senor  provide  high  spatial  resolutions.  Using  this  high  performance  thermometer,  excellent  characterization  of  the  real¬ 
time  temperature  profile  within  the  flow  of  water  turbulence  has  been  realized. 

Keywords:  Fiber-optic  sensors,  Fabry-Perot  interferometer,  thermometer,  high-speed  detection. 


1.  INTRODUCTION 

Temperature  measurement  is  an  important  aspect  of  oceanography,  which  functions  as  a  key  step  toward  field 
application  in  both  civilian  and  military  environment  [1,  2].  Currents,  mixing  events,  large-scale  circulation,  air-sea 
exchange,  and  certainly  global  climate  change  are  a  few  obvious  examples.  Temperature  fluctuations,  especially  on  the 
small  scale,  also  affects  sensor  and  system  performance.  For  example,  diver  visibility,  search  and  rescue  at  sea,  mine 
detection  all  rely  on  the  underwater  imaging  capability.  It  has  been  shown  recently,  that  not  only  the  turbidity  affects 
image  quality,  but  also  optical  turbulence.  As  we  know,  turbidity  is  caused  by  particle  absorption  but  mostly  scattering  in 
the  water.  It  thus  may  be  inferred  that  this  issue  does  not  exist  in  clean  waters  like  clean  lakes  and  rivers.  But  in  reality, 
another  factor,  i.e.,  the  optical  turbulence,  can  take  over  and  play  the  dominant  role  in  determining  the  imaging  quality  in 
these  clear  waters  [3],  and  can  also  affect  acoustical  signal  propagation  [4]. 

While  salinity  variations  could  sometimes  lead  to  severe  turbulence  [5],  temperature  gradient  is  the  primary 
contributing  element  that  triggers  the  underwater  optical  turbulence  by  modifying  the  index  of  refraction  within  the 
microstructures.  To  better  understand  the  process,  and  possibly  to  facilitate  rectification  of  the  image  blurred  by  the 
turbulence,  a  full  characterization  of  the  temperature  distribution  at  the  highest  sampling  rate  and  at  the  smallest  spatial 
scale  is  critical.  Specifically,  it’s  reported  that  a  sensor  with  singe-pole  time  constant  less  than  and  spatial  resolution 
better  than  is  required  to  estimate  the  rate  of  dissipation  of  temperature  variance  in  turbulent  heat  flux  model  [6].  To 
tackle  the  issue  regarding  the  fast  and  high-resolution  measurement  of  localized  temperature,  several  efforts  have  been 
made  by  resort  to  a  tiny  thermistor  or  thermocouple  [6,  7].  However,  interferences  from  the  environment,  such  as  strong 
electromagnetic  field  from  lightening,  erroneous  signals  or  corrosion  from  strong  salinity,  etc.,  prevent  the 
thermocouples  from  being  reliable  sensors.  Furthermore,  the  long-distance  transmission  of  the  weak  electronic  signals  is 
a  big  challenge. 

In  this  paper,  we  will  study  an  alternative  high-resolution  and  fast-response  fiber-optic  temperature  sensor  we 
recently  developed  [8].  It  consists  of  a  tiny  silicon  cylinder  assembled  on  the  endface  of  a  single-mode  fiber  (SMF).  The 
silicon  cylinder  serves  as  a  Fabry-Perot  interferometer  (FPI)  which  sends  back  an  interferometric  spectrum  that  is 
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sensitive  to  the  silicon  temperature.  Due  to  the  large  thermo-optic  coefficient  of  silicon,  the  investigated  sensor  has  a 
high  temperature  sensitivity.  To  increase  the  temperature  resolution  by  reducing  noise,  an  average  wavelength  tracking 
previously  developed  by  us  will  be  applied  [8].  The  small  size  of  the  sensing  silicon,  no  longer  than  200  pm  and  no 
thicker  than  100  pm,  endows  the  sensor  with  miniature  footprint  which  is  a  big  advantage  for  surveillance  of  temperature 
variation  at  a  tiny  spot.  With  such  small  size,  along  with  the  high  thermal  diffusivity  of  silicon,  the  fiber  sensor  features 
swift  response  to  follow  the  fast  temperature  variation  within  the  turbulence.  Furthermore,  the  extremely  small  optical 
attenuation  (0.1  dB/km)  within  the  SMF  is  superior  when  it  comes  to  remote  sensing  in  far  and  deep  ocean.  As  an  initial 
test,  the  real-time  temperature  structure  within  the  water  microstructure  generated  by  two  methods  will  be  demonstrated. 
Comparison  in  performance  of  the  investigated  sensor  with  the  mostly  used  commercial  high-speed  thermistor  FP07  (GE 
Thermometries)  will  also  be  provided. 


2.  PRINCIPLE  AND  THEORY 


Figure  1.  (a)  Schematic  structure  of  the  sensor  head  and  (b)  illustration  of  the  principle  of  operation. 


Structure  of  the  sensor  head  is  shown  in  Fig.  1(a).  According  to  the  procedure  described  in  detail  in  our  previous  report 
[8],  the  cylindrical  silicon  can  be  mounted  on  the  cleaved  endface  of  a  SMF.  The  silicon  cylinder  works  as  a  Fabry-Perot 
interferometer  (FPI),  as  shown  in  Fig.  1(b).  Two  reflections  at  the  fiber-silicon  and  silicon- surrounding  interfaces 
combine  and  send  back  an  interferometric  spectrum  which  features  a  series  of  dip  and  peak  wavelengths.  If  refractive 
index  of  the  silicon  material  is  defined  as  n  and  length  of  the  cylinder  is  denoted  as  L,  then  the  optical  path  difference 
(OPT)  between  these  two  reflected  lights  will  determine  the  Nth  peak  wavelength,  AN,  following  this  formula  below 

(N+Y7)XN=2nL.  (1) 


The  inherent  properties,  i.e.,  thermo-optic  effect  and  thermal  expansion,  give  birth  to  the  temperature  dependence  of 
both  n  and  L  in  Eq.  (1).  Consequently,  the  peak  wavelength  AN  is  a  function  of  temperature,  which  is  explored  to 
monitor  the  environmental  temperature  variation.  Temperature  sensitivity,  k ,  of  the  peak  wavelength  derived  from  Eq. 
(1)  can  be  expressed  as 


i  i  ( 1  1 

k  — - —  An - b - 

dT  \ndT  L  dT ) 


(2) 


Right-hand  side  (RHS)  of  Eq.  (2)  includes  two  parts,  one  is  the  contribution  of  thermo-optic  effects,  kT0  = 
the  other  part  is  contribution  of  thermal  expansion,  kTE  =  As  an  estimation,  kT0  =  68Apm/°C  and  kTE  = 

4pm/°C  are  obtained  using  the  thermo-optic  coefficient  (TOC)  of  1.5  x  10-4  RIU /°C  and  thermal  expansion 
coefficient  (TEC)  of  2.55  x  10-6  m/(m  ■  °C)  at  temperature  of  25  °C,  evaluated  around  1550  nm  and  n  =  3.4  for 
silicon.  It’s  clear  that  the  total  sensitivity  of  k  =  kT0  +  kTE  =  72Apm/°C  is  mainly  attributed  to  the  thermo-optic 
effects. 

To  predict  the  speed  of  this  sensor  in  response  to  temperature  change,  a  heat  transfer  analysis  should  be  performed. 
As  our  previous  work  has  provided  theoretical  analysis  on  this  issue  [8],  a  brief  reproduction  of  the  theory  is  extracted 
here  to  facilitate  interpretation  and  evaluation  of  the  sensor  response  time  with  extended  detailed  analysis  on  the 
influence  of  sensor  size  on  the  response  time.  Due  to  the  high  thermal  conductivity  and  small  thermal  mass,  whenever 
there  is  a  thermal  exchange  between  the  sensor  head  and  surrounding,  the  heat  energy  dissipates  so  fast  within  the  silicon 
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cylinder  that  the  temperature  gradient  is  ignorable.  In  such  a  situation,  the  formula  for  heat  transfer  analysis  is  reduced  to 
the  following  lamped  parameter  model  [9] 

hAs(Tf  ~T)dt  =  pC  VdT ,  (3) 


where  T ^  and  T  are,  respectively,  the  steady-state  and  transient  temperature  of  silicon  cylinder;  As  is  the  surface  area 
through  which  heat  exchanges;  ps,  Cs,  and  Vs  are,  respectively,  mass  density,  heat  capacity  and  volume  of  the  silicon 
FPI;  h  is  the  heat  transfer  coefficient.  Solution  to  Eq.  (3)  is  given  as 

fhA 

T(t)  =  Tf  +(Tn-  Tf ) exp 


(4) 


v  P.C.V,  , 

where  T0  is  the  initial  temperature  of  the  silicon  sensor  head.  Thus,  combination  of  Eqs.  (2)  and  (4)  produces  the 
following  temperature-dependent  peak  wavelength 

(  hA 

A(t)  =  T0  +  kTx  +  & (f0  - 7^ ) exp - 


(5) 


v  pscys  j 

where  A0  is  the  intercept  wavelength.  Obviously,  the  RHS  of  Eq.  (5)  includes  a  transient  item,  from  which  the  response 
time,  defined  as  the  time  needed  when  the  signal  rises  to  1/e  of  the  total  change,  is  extracted  as  follows 

T=P'C'V‘  ■  (6) 

hAs 

To  obtain  the  value  of  response  time  t,  it’s  necessary  to  first  get  the  value  of  heat  transfer  coefficient  h.  Using  the 
Nusselt  number  Nu  and  heat  conductivity  Kw  of  water,  h  can  be  calculated  as  [9] 


(7) 


where  Lc  =  VS/AS  is  the  characteristic  length.  Substituting  Eq.  (7)  into  Eq.  (6)  and  taking  As  =  nDL  +  nD2 /4  and 
Vs  =  nD2L/A ,  where  D  and  L  are,  respectively,  diameter  and  length  of  the  silicon  cylinder,  give  rise  to 

p,c,  l2d 2 


NuKw(4L  +  D)2 

Eq.  (8)  gives  an  explicit  dependence  of  response  time  on  the  size  of  the  silicon  cylinder. 


(8) 


Figure  2.  Calculated  response  time  as  a  function  of  (a)  diameter  and  (b)  length  of  the  cylindrical  silicon  FPI. 


Taking  Nu  =  1.2  for  the  heat  transfer  on  the  water/silicon  surface  [10],  Kw  =  0.6  W/(m  ■  K),  Cs  =  712  J/(kg  ■  K), 
and  ps  =  2.329  x  103  kg/m3,  the  calculated  response  time  versus  diameter  and  length  of  silicon  cylinder  is  shown  in 
Fig.  2.  It  can  be  inferred  from  Fig.  2(a)  that  by  reducing  diameter  of  the  silicon  cylinder  the  response  time  decreases. 
Also,  by  reducing  length  of  the  silicon  cylinder  the  response  time  decreases,  as  illustrated  in  Fig.  2(b).  If  diameter  and 
length  of  the  silicon  pillar  are,  respectively,  50  pm  and  100  jam,  then  a  response  time  as  short  as  -0.3  ms  is  evaluated. 
The  reduced  size  (diameter  and  length)  of  the  sensing  silicon  facilitates  the  heat  exchange  between  the  sensor  head  and 
environment,  therefore,  it  is  more  favorable  if  reduced  sensing  mass  is  applied.  However,  the  reduced  size  of  the  sensor 
head  will  extremely  increase  the  difficulty  of  the  assembling  procedure  depicted  in  detail  in  [8].  Furthermore,  the 
divergent  outgoing  rays  from  the  single-mode  fiber  may  expand  out  of  the  silicon  pillar  if  its  diameter  is  too  small, 
leading  to  a  degraded  visibility  and  thus  limited  temperature  resolution.  Lastly,  we  have  demonstrated  in  [8]  that  an 
average  wavelength  tracking  method  taking  advantage  of  the  dense  fringes  will  help  reduce  noise,  but  a  reduced  length 
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of  the  silicon  cylinder  is  a  contradiction  to  this  advantage,  because  a  reduced  length  results  in  a  larger  free  spectral  range 
and  consequently  fewer  peaks  used  for  calculating  the  average  wavelength.  Therefore,  a  balance  between  response  time 
and  resolution  should  be  considered  while  choosing  the  size  of  the  sensing  element. 


3.  EXPERIMENTATION 
3.1  Temperature  resolution  and  response  time 


Computer  High-speed 

Spectrometer 


>er Probe 

Water 


Figure  3.  Experimental  setup 

Experimental  demonstration  of  the  sensor  performance  has  been  carried  out  using  the  setup  exhibited  in  Fig.  3.  The 
broadband  source  provides  the  wide-band  spectrum  for  average  wavelength  tracking.  The  high-speed  spectrometer  (I- 
MON  256,  Ibsen  Photonics)  with  a  maximum  sampling  frequency  of  6  kHz  guarantees  the  swift  capture  of  spectrum 
variation  over  time. 
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Figure  4.  (a)  A  typical  reflection  spectrum  obtained  from  the  sensor,  (b)  Relative  wavelength  shift  as  a  function 
of  temperature,  (c)  Wavelength  noise  versus  time,  (d)  Temporal  wavelength  while  dipping  the  sensor  head  from 

air  (-20  °C)  into  hot  water  (-58  °C). 
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One  typical  frame  of  the  spectrum  is  shown  in  Fig.  4(a),  from  which  dozens  of  fringe  peaks  are  inspected,  owing  to 
the  small  free  spectrum  range  of  only  1.726  nm.  To  make  full  use  of  the  large  amount  of  peaks,  we  have  proposed  and 
investigated  an  average  wavelength  tracking  method  to  reduce  the  noise  [8],  which  is  also  applied  in  this  paper.  In  the 
experimental  results  shown  in  Fig.  4,  we  used  22  peaks  that  were  centered  around  1545  nm.  The  temperature  sensitivity 
was  first  measured  and  the  results  are  shown  in  Fig.  4(b).  An  excellent  linear  response  was  observed  within  the 
temperature  range  of  20-100  °C,  indicating  a  temperature  sensitivity  of  84  pm/°C.  By  immersing  the  sensor  into  the  water 
stabilized  in  the  room  temperature  (around  22  °C)  and  recording  the  spectrum  for  one  second  with  a  sampling  frequency 
of  6  kHz  and  integration  time  of  100  jus,  the  wavelength  noise  was  characterized  and  the  results  are  shown  in  Fig.  4(c). 
The  standard  deviation  of  -0.0504  pm  of  the  wavelength  variation  during  this  period  is  defined  as  the  noise  level, 
leading  to  a  temperature  resolution  of  6><10"4  °C.  The  high  resolution  guarantees  the  visualization  of  subtle  variation  in 
the  local  water.  To  test  the  response  time  of  the  proposed  sensor,  the  sensor  head  was  fast  dipped  into  a  cup  of  hot  water. 
The  temporal  variation  in  wavelength  is  shown  in  Fig.  4(d),  from  which  it  can  be  seen  that  the  fast  response  of  the  sensor 
started  to  show  up  before  it  hit  the  water  at  the  time  of  around  3.5  ms,  which  was  due  to  the  warmed  up  air  and  water 
vapor.  After  that  a  sharp  increase  in  the  wavelength  was  observed  upon  immersing  in  the  hot  water.  From  the  curve  in 
Fig.  4(d),  a  response  time  of  0.5  ms  is  inferred,  suggesting  a  sampling  frequency  as  high  as  2  kHz  can  be  achieved. 

3.2  Flow  sensing  in  water  turbulence 

To  demonstrate  the  performance  of  the  investigated  sensor  in  profiling  the  temperature  variation  within  water  turbulence, 
two  experiments  were  carried  out.  In  the  first  experiment,  the  turbulence  was  created  by  putting  an  icebag  on  the  surface 
of  water.  The  water  cooled  down  by  the  icebag  went  down  quickly  and  in  the  meantime  warm  water  rose,  thus  creating  a 
turbulence  and  leading  to  the  variation  in  temperature  at  the  fixed  point  where  the  fiber  sensor  was  placed.  The  whole 
process  was  perfectly  monitored  by  our  fiber  thermometer,  as  shown  by  Fig.  5(a). 

In  the  second  experiment,  the  turbulence  was  generated  by  pouring  some  hot  water  into  a  bucket  of  cool  water.  This 
time,  the  fiber  sensor  was  compared  with  a  commercial  fast-response  thermistor  (FP07)  which  served  both  as  a  reference 
and  as  a  calibration.  The  two  sensors  were  deployed  in  the  vicinity  of  each  other.  Comparison  between  these  two  sensors 
is  shown  in  Fig.  5(b).  Apparently,  the  fiber  sensor  followed  the  FP07  thermistor  very  well  but  tracked  much  more  details 
of  the  swift  temperature  alteration.  These  useful  temperature  profile  will  provide  a  powerful  tool  for  the  analysis  of 
influence  of  underwater  turbulence  on  imaging  in  the  ocean.  Data  were  also  collected  in  a  simulated  turbulence 
environment,  utilizing  a  Rayleigh-Benard  convective  tank.  More  details  can  be  found  in  a  companion  paper  by  Matt  et 

al.  [11]. 


Figure  5.  Detection  of  water  turbulence  generated  by  (a)  putting  an  ice  bag  on  the  surface  of  water  and  (b)  pouring  some  hot 
water  into  a  bucket  of  cool  water.  In  (b),  response  of  the  commercial  high-speed  thermistor  FP07  is  provided  as  calibration 

and  reference. 

4.  CONCLUSION 

In  summary,  this  work  describes  a  novel  fiber-optic  thermometer  that  features  high  temperature  resolution  and  fast 
response  for  the  oceanographic  applications.  Theoretical  investigation  suggests  that  an  increase  in  sampling  rate  can  be 
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obtained  by  reducing  the  size  (diameter  or  length)  of  the  silicon  cylinder.  For  example,  a  50-jam-diameter  and  100-jam- 
length  silicon  pillar  is  theoretically  predicted  to  possess  a  response  time  as  low  as  0.3  ms.  In  experiments,  a  sensor  with 
80-jam-diameter  and  200-pm-length  silicon  pillar  was  capable  of  resolving  a  temperature  variation  as  low  as  6><10"4  °C 
and  sampling  at  a  frequency  as  high  as  2  kHz.  Using  this  sensor,  temperature  variation  within  turbulence  micro  structure 
has  been  successfully  tracked.  The  water  turbulence  was  initiated  by  two  ways,  one  is  floating  an  icebag  on  the  water, 
and  the  other  is  pouring  some  hot  water  into  a  bucket  of  cool  water.  Outperformance  of  the  proposed  fiber  sensor  over  a 
commercial  high-speed  thermistor  FP07  in  characterizing  the  temperature  fluctuation  within  the  turbulence  manifests 
this  powerful  tool,  which  exhibits  a  great  promise  for  applications  in  oceanography. 
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